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Abstract. The objective of this paper is to develop a model and a mini- 
mization method to provide flight path optimums reducing aircraft noise in 
the vicinity of airports. Optimization algorithm has solved a complex optimal 
control problem, and generates flight paths minimizing aircraft noise levels. 
Operational and safety constraints have been considered and their limits sat- 
isfied. Results are here presented and discussed. 
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Nomenclature 



c.„ 


slope of lift coefEcient curve 


(^Xo 


drag coefficient 


T 


thrust, N 


L 


lift, N 


D 


drag, N 


9 


9.8 m/s^ 


X 


horizontal distance, m 


y 


lateral distance, m 


h 


aircraft height, m 


k 


induced drag coefEcient 


t 


time, s 


M 


roll angle, rad 


a 


angle of attack, rad 


7 


flight path angle, rad 


X 


yaw angle, rad 



5x 


throttle setting 


h 


induced drag parameter 


V 


speed of aircraft, m/s 


s 


wing area, m^ 


c 


speed of sound, m/s 


M 


mach number, V/c 


P 


air density, kg/m^ 


d 


nozzle diameter, m 


s 


area of coaxial engine nozzle, m 


V 


speed of gas, m/s 


w 


density exponent 


T 


temperature, K 


e 


directivity angle, rad 


R 


source-to-observer distance, m 


m 


aircraft mass, kg 



X 



1. Introduction 

Since the introduction of jet aircraft in the 1960, aircraft noise produced in the 
vicinity of airports has represented a serious social and environmental issue. It is 
a continuing source of annoyance in nearby communities. The importance of that 
problem has been highlighted by the increased public concern for environmental 
issues. To deal with this problem, aircraft manufacturers and public establishments 
are engaged in research on technical and theoretical approaches for noise reduction 
concepts that should be appHed to new aircraft. The ability to assess noise exposure 
accurately is an increasingly important factor in the design and implementation of 
any airport improvements. 
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Aircraft are complex noise sources and the emitted intensities vary with the type 
of aircraft, in particular, with the type of engines and with the implemented flight 
procedures. The noise contour assessment due to the variety of flight route schemes 
and predicted procedures is also complex. A set of data must be used which includes 
noise data, flight path parameters and their features, and environmental conditions 
affecting outdoor sound propagation. Three development initiatives are available 
for the reduction of aircraft noise: (1) innovative passive technologies required by 
the industry for developing environmentally compatible and economically viable air- 
craft, (2) advanced active technologies such as computational aeroacoustics, active 
control, advanced propagation and prediction methods, (3) reliable trajectory and 
procedures optimization which can be used to determine optimal landing approach 
for any arbitrary aircraft at any given airport. The last action will be particularly 
emphasized in the next sections. 

A number of calculation programs of aircraft noise impact have been developed 
over the last 30 years. They have been widely used by aircraft manufacturers and 
airport authorities. Their rehability and results efficiency to assess the real impact 
of aircraft noise have not been proved conclusively. They are complex, very slow 
and can not be planned for on-Hne and on-site use. That is the reason why, the 
model described in this paper, generating optimal trajectories minimizing noise, is 
considered as a promising scientiflc plan. 

Several optimization codes for NLP exist in the literature. After a large number of 
test and comparisons, we choose KNITRO [7] which known for its performances and 
robustness, this software is efflcient to solve general nonlinear programming. We 
will explain in the following sections how the considered optimal control problem 
is discretized and solved. Numerical results have been analyzed and their relia- 
bility and flexibility have been proved. We have demonstrated the effectiveness 
of computation and its application to aircraft noise reduction. The objective of 
that alternative research is to develop high payoff models to enable a safe, and 
environmentally compatible and economical aircraft. We should make large proflts 
in terms of noise abatement in comparison with the expected noise control sys- 
tems in progress. These systems, which are not in an advanced step, in particular 
at low frequencies, are still ineffective or impractical. Actually, the low-frequency 
broadband generated by the engines represents a signiflcant source of environmen- 
tal noise. Their radiation during flight operations is extremely difficult to attenuate 
using the mentioned systems and is capable of propagating over long distance |21] . 

Details of trajectory and aircraft noise models, and optimal control problems are 
presented in section 2 and 3 while the last section is devoted to numerical experi- 
ments. 



2. Optimal Control Problem 

2.1. Equation of Motion. In general, the system of differential equations com- 
monly employed in aircraft trajectory analysis is the following six-dimension system 
derived at the center of mass of aircraft : 
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where T = r(/i, F, 4), D = D(h, V, a) and L = L{h, V, a). 

These equations embody the assumptions of a constant weight, symmetric flight 

and constant gravitational attraction [U [6] . 

Figure (1) shows the forces acting on an aircraft at its center of gravity during 
an approach. 
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Figure 1. Aircraft forces during phase of descent 



Those equations could be applied to conventional aircraft of all sizes. The most 
dominant aerodynamics affecting results are the lift L and drag D, deflned as follows 

El 

D = ysv^ [a,, + hciy] 

The thrust model chosen, by Matthingly [l^, depends explicitly on the aircraft 
speed, the geometric aircraft height and the throttle setting 

Po V 2 

To is full thrust, po is atmospheric density at the ground (= 1.225 kg/rrt') and p is 
atmospheric density at the height h{p — p^{l — 22.6 x 10~^ft.)'*-^^). 



4 OPERATIONAL AIRCRAFT PARAMETERS REDUCING NOISE EMISSION 

The previous system of equations (|2.ip can be written in the following generic 
form : 

X{t) = f{X{t),U{t)) 

where : 

X : [to,tf] -^ M^ 

t — > X{t) = [F(i),7(t),x(i),a;(i),y(i), /i(i)] are the state variables, 

U: [to,tf] -^ m"" 

t — > C/(i) = [Q;(t), (5a; (i),/i(i)] are the control variables, 

and io and t/ are the initial and final times. 



2.2. Constraints. Search for optimal trajectories minimizing noise must be done 
in a realistic flight domain. Indeed, operational procedures are performed with re- 
spect to parameter limits related to the safety of flight and the operational modes 
of the aircraft. 

• The throttle stays in some interval 

• The speed is bounded 

l.Ws, <V< Vmax 

where Vgg the stall velocity, the limited velocity at which the aircraft can 
produce enough lift to balance the aircraft weight. Vgg and Vmax depend 
on the type of aircraft. 

• The flight path angle providing a measure of the angle of the velocity to 
the inertial horizontal axis, is bounded 

^niin _ T — '~imax 

• The angle of attack is bounded 

O^niin Jl ^ 2l ^max 

• The yaw angle and roll angle stay in some prescribed interval 



Those inequality constraints could be formulated as: 

a< C{X{t),U{t)) < b 

where 

C: M^xM^ — > M^ 

^^> {X{t),Uit)) -^ C{X{t),Uit)) = [^it),V{t),x{t),a{t),SAt),Kt)] 

a and b are two constant vectors of IR^. 
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2.3. Cost function. Models and methods used to assess environmental noise prob- 
lems must be based on the noise exposure indices used by relevant international 
noise control regulations and standards (ICAO O [TO], Lambert and Vallet [9]). 
As described by Zaporozhest and Tokarev p9], these indices vary greatly one from 
another both in their structure, and in the basic approaches used in their definitions. 

The cost function to be minimized may be chosen as any usual aircraft noise in- 
dex, which describes the effective noise level of the aircraft noise event [IHl [30] , like 
SEL (Sound Exposure Level), the EPNL (Efi^ective Perceived Noise Levels) or the 
Leq,At (Equivalent noise level), ... It is well known that the magnitude of Leg, At 
correlates well with the effects of noise on human activity, in particular, with the 
percentage of highly noise-annoyed people living in regions of significant aircraft 
noise impact. This criterion is commonly used, as basis, for the regulatory basis in 
many countries. Based on comparison of noise exposure indices and a comparison 
of the methodologies used to calculate the aircraft noise exposure, it can be con- 
cluded that the general form of the most used and accepted noise exposure index 
is Leq.At that we have chosen as an index (|1HI15)). L^q^AT is expressed by: 

(2) ie,,AT =101og^y 'lO°l^-(*)di 

where io initial time, i/ final time, AT ^ tf — to and Lp{t) is the overall sound 
pressure level (in decibels). 

We will define in the next subsection the analytic method to compute the noise 
level at any reception point. 

Calculation method for aircraft noise levels 

The aircraft noise levels Lp at a receiver is obtained by the following formula 
based on works [TBllSTI: 



(3) Lp = Lref - 20 logio R + ^atm + Aground + Ay + Ap + A/ 

where Lref is the sound level at the source, 20 logio ^ is a correction due to geo- 
metric divergence, Aatm is the attenuation due to atmospheric absorption of sound. 
The other terms Aground, Ay, A^, A/ correspond respectively to the ground ef- 
fects, correction for the Doppler, correction for duration emission and correction 
for the frequency. 

In this paper, we have used a semi-empirical model to predict noise generated 
by conventional- velocity-profile jets exhausting from coaxial nozzles predicting the 
aircraft noise levels represented by the jet noise [22] which corresponds to the main 
predominated noisy source. It is known that jet noise consists of three principal 
components. They are the turbulent mixing noise, the broadband shock associated 
noise and the screech tones. At the present time, this first approximation have been 
used herein. It seems to be correct in that step of research because the complexity 
of the problem. Many studies have agreed with this model and full-scale experi- 
mental data even at high jet velocities in the region near the jet axis. Numerical 
simulation of jet noise generation is not straightforward undertaking. Norum and 
Brown j^, Tam and Auriault [23l[2lll2S] had earher discussed some of the major 
computational difficulties anticipated in such effort. At the present time, there are 
reliable to jet noise prediction. However, there is no known way to predict tone 
intensity and directivity; even if it is entirely empirical. This is not surprising for 
the tone intensity which is determined by the nonlinearities of the feedback loop. 
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Obviously, to complete this study we will need to integrate other noise source mod- 
els in particular aerodynamics. 



Although the numerous aspects of the mechanisms of noise generation by coax- 
ial jets are not fully understood, the necessity to predict jet noise has led to the 
development of empirical procedures and methods. During the descent phase, the 
jet aircraft noise as well as propeller aircraft noise is approximately omni-directional 
and the noise emission is decreasing with decreasing speed when assuming that the 
power setting is constant. The jet noise results from the turbulence created by 
the jet mixing with the surrounding air. Jet mixing noise caused by subsonic jets 
is broadband in nature (its frequency range is without having specific tone com- 
ponent) and is centered at low frequencies. Subsonic jets have additional shock 
structure- related noise components that generally occur at a higher frequencies. 
The prediction of jet noise is extremely complex. The used methods in system 
analysis and in the engine design usually employed simpler or semi-empirical pre- 
diction techniques. By replacing the predicted jet noise level [22j in |[3l), we obtain 
the following expression 
(4) 
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141 + 101ogio(^ 
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v\ 

V2 
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T2 
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-201ogioi?- 151ogio(CD(Afc,0)) - 101ogio(l -Mcos^), 

speed of jet gas at inner contours 

speed of jet gas at outer contours 

area of coaxial engine nozzle at inner contours 

area of coaxial engine nozzle at outer contours 

temperature at inner contours 

temperature at outer contours 

atmospheric density at inner contours 

International Standard Atmosphere density 

International Standard Atmosphere for speed of sound. 



The effective speed T4 is defined by: 

V, = vi[l-{V/vi)cos{aT)f'''. 

The angle of attack ax, the upstream axis of the jet relative to the direction of 
aircraft motion has been neglected in this study. 
The distance source to the observer is 

R={X- Xobsf + {y - Vobsf + h"^ 

where {xobs,yobs) are the coordinates of the observer. Ay is expressed by 
Ay = -lb\og{CD{Mc,6)) - 101og(l - Mcose) 
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where Cd{Mc.,6) indicate the Doppler convection factor: 

CDiMc,e) = [(1 + Afc cos 61)2 + 0.04M2], 
the Mach number of convection is: 

Mc = 0.62(wi - V)/c 

and w — ; , ,„ ^ — 1. 

o.6 + (v;/c)3-5 

The validity of this improved prediction model is established by fairly extensive 
comparisons with model-scale static data ^26j. InsufRcient appropriate simulated- 
flight data are available in the open literature, so verification of flight effects during 
aircraft descent has to be established. Analysis by Stone and al. (22] has shown 
that measured data are used to calibrate the behavior of the used jet noise model 
and its implications from theory. Nowadays, the above formulation is being consid- 
ered realistic compared to others described in applied and fundamental literature 
deahng with jet noise of aircraft during descent operations. 

Taking account the formulas ^ and ^ , we obtain our cost function in the following 
integral function form 



{Xit),U{t)) -^ J{X{t),U{t))^ i{X{t),U{t))dt 



J is the criterion for optimizing the noise level at the reception point. It doesn't 
depend of Xj 7 and U. 

2.4. Optimal Control Problem. Finding an optimal trajectory, in term of min- 
imizing noise emission during a descent, can be mathematically stated as an ODE 
optimal control problem. We opted different notations z = X,u = U and to — 0: 

min J{z,u) = / £{z{t),u{t))dt 



(OCP) 



Jo 

z{t) ^ f{z{t),u{t)),yt e[0,tf] 

Z/l(0) = Cl, Zl^{tf) = C2 

a < C{z{t),u{t)) < b 

where J : 5?"+™ -^ IRJ : 1R''+"' -^ R" and C : K"+" -^ If correspond 
respectively to the cost function, the dynamic of the problem, and the constraints 
defined in the previous section. The initial and final values for the sate variables 
((a;(0),y(0),/i(0),y(0)) and {x{tf),y{tf),h{tf)) are fixed, m, := \h\ + {hi is the 
total number of fixed limit values of the state variables. 

To solve our problem (find the optimal control u{t) and the corresponding opti- 
mal state z(t)), we discretize the control and the state with identical grid and 
transcribe optimal control problem into nonlinear problem with constraints. The 
next section may be helpful in telling how the problem could be solved. They 
will present theoretical consideration and computational process that yield to fiight 
paths minimizing noise levels at a given receiver. 
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3. Discrete Optimal Control Problem 

To solve {OCP) different methods and approaches can be used [Sjl^Sj. In this 
paper, we use a direct optimal control technique : we first discretize {OCP) and 
then solve the resulting nonlinear programming problem. 

3.1. Discretization. We use an equidistant discretization of the time interval as 

= io < ••■ < ^Af = tf 
where : 



tk=tQ + kh, fc = 0, ..., N and h 



N ■ 



Then we consider that u{.) is parameterized as a piecewise constant function : 

u(t) := Uk for t e [tk-i,tk[ 



and use a Runge-Kutta scheme (Heun) to discretize the dynamic : 

s 

Zk+1 = Zk + \i^bif(zki,Uk) 
1=1 

s 

Zki = Zfc + h ^ aijf{zkj , Uk) 
k = 0,...,N ~1, i = l,...,s. 
The new discrete objective function is stated as : 

N 

J = y^j{zk,uk). 

k=0 

The continuous optimal control problem {OCP) is replaced by the following dis- 
cretized control problem : 

N 

miny^^(zfc,Mfc) 

fe=0 



{NLP) 



Zk+i = Zk +h^6j/(zfej,Ufe), A: = 0, . ..,iV- 1 



4=1 



Zki = Zk +h^ay/(zfcj,Ufe), fc = 0,. ..,7V- 1 

ZOi^ = Ci , ZNj^ = C2 

. a<C{zk,Uk) <6, fc = 0,...,iV 



To solve (NLP) we developped an AMPL [2] model and used a robust interior point 
algorithm KNITRO [7]. We choose this NLP solver after numerous comparisons 
with some other standard solvers available on the NEOS (Server for Optimization) 
platform. 

4. Numerical results 

For different cases and configurations, we consider an aircraft approach with an 
initial condition {xq — 0;yQ = 0; Hq — 3500 to) a final condition (x/ = 60000 to,; yf = 
5000 m; hf — 500 to) and for a fixed tf — 10 min and a discretization parameter 
TV = 100 or TV = 200. 
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We first consider the simplest configuration of one single observer and no additional 
constraint. 

4.1. One fixed observer. For various positions {xobs,yobs) of an observer on the 
ground (near the aircraft trajectory) we calculate the optimal noise level J (corre- 
sponding to our optimal trajectory Tr) and the noise level Ji corresponding to the 
trajectory Tri that minimizes the "fuel consumption" (re minimizing) the simple 
following model of consumption P : 



C0{h,V,5,)^ / CsRT(t)dt 
Jo 
where Csr is supposed constant. 



The following table summarizes the obtained results. 



y^obs 1 yobs) 


J{dB) 


max{f.e, o.e) 


CPU (s) 


Ji{dB) 


Ji - J{dB) 


%C0 


(0,0) 


45.92 


7.92e-07 


10 


47.03 


1.1 


36% 


(0,2500) 


44.95 


8.76e-07 


8.4 


46.32 


1.4 


36% 


(0,5000) 


43.27 


3.44e - 07 


9.8 


44.93 


1.7 


36% 


(20000,0) 


48.97 


4.66e - 07 


10.4 


51.04 


2.10 


33% 


(20000, 2500) 


49.58 


2.81e-08 


10.6 


51.55 


2 


26% 


(20000, 5000) 


47.18 


8.36e-07 


9.8 


49.70 


2.5 


36% 


(40000, 0) 


47.59 


6.73e-07 


14.6 


49.52 


2 


26% 


(40000, 2500) 


50.76 


4.21e-07 


6.9 


52.87 


2.11 


26% 


(40000, 5000) 


49.74 


6.92e-07 


11.9 


51.72 


2 


24% 


(60000,0) 


42.85 


5.29e-07 


7.2 


45.00 


2.15 


34% 


(60000, 2500) 


45.04 


6.90e-07 


6.7 


48.014 


3 


36% 


(60000, 5000) 


48.84 


7.10e-07 


6 


54.18 


5.34 


38% 



Table 1. Noise minimization 



For each case, the algorithm (KNITR0[7j) found a solution with a very high accu- 
racy. The computation of Tri have been done only one time; it needs 7s with an 
accuracy of max{f.e, o.e) = 4.24e — 07. 

The third column of Table 1 measure the maximum of feasibility error and op- 
timality error, the fourth one gives an idea on the computation effort (namely the 
CPU time). The two last columns correspond to the noise reduction and % of 
exceeded consumption : %C0 = {CO{Tr) - CO{Tri))/CO{Tr). 
Our trajectory that minimizes the noise consume about 31% more than the trajec- 
tory minimizing the consumption. 



The following figure showing the solution trajectory Tr, where the fixed observer 
presents a certain area near the airport : 



10 
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5000 



Figure 2. Trajectory in 3D 

The following figures present the state and control variables of the optimal tra- 
jectory Tr. We remark that the optimal variables ft-, V,x,a and /i present some 
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Figure 3. Solutions of {NLP) 



large constant stages, while 7 and 5^ are bang-bang. 

4.1.1. One fixed observer with an additional consumption constraint. Table (1) 
shows that the optimal trajectory Tr consumes about 31% more that Tri. This 
fact makes of interest some additional constraint on the consumption. 
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We define a new problem : 
(5) {0CP)2 



min Leq,At 
i(t)-/(z(i),«(i)) 

Z/i(0) = Cl, Zl^{tf) = C2 

CO{z{t),u{t)) <= l.lCO{Tri) 



This problem can be solved using the same techniques (discretization,...) and the 
same configurations. We obtain the following results. 



{Xobs,yobs) 


J{dB) 


max{f.e, o.e) 


MdB) 


CPU(s) 


(0,0) 


46.32 


9.94e-07 


47.03 


9.6 


(0,2500) 


45.44 


9.17e-07 


46.32 


13.2 


(0,5000) 


43.94 


9.49e-07 


44.93 


13.1 


(20000,0) 


49.79 


7.25e-07 


51.04 


18.6 


(20000, 2500) 


50.32 


8.75e-07 


51.55 


8.5 


(20000, 5000) 


48.33 


9.37e-07 


49.70 


20 


(40000, 0) 


48.16 


7.54e-07 


49.52 


33.3 


(40000, 2500) 


51.43 


6.55e-07 


52.87 


9.3 


(40000, 5000) 


50.31 


6.90e-07 


51.72 


20.4 


(60000,0) 


43.81 


9.74e-07 


45.00 


34.8 


(60000, 2500) 


46.40 


7.69e-07 


48.014 


27.9 


(60000, 5000) 


51.32 


9.30e-07 


54.18 


17.8 



Table 2. Noise minimization 



The following figures present the state and control variables and the solution tra- 
jectory Tr : 




Figure 4. Trajectory in 3D 
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Figure 5. Solutions of {0CP)2 

We obtain approximately the same characteristic for the trajectory while the 
noise reduction is still significative. 



4.2. Several observers fixed on the ground. We can generalize the processus 
of minimization for several observers. In this case, we minimize the maximum of 
noise corresponding to several observers. 
The problem to solve is written as follows : 



(6) 



[OCP), 



mini? 

■d >^ Jobs, 

Z/i(0) = Ci, Zi^{tf) = C2 



where Jobsj are the noise levels corresponding to j fixed observers. 

Once again, we use a direct method to solve the problem (0CP)3 with the same 
modeling language and software. 

We choose five observers : (o5si(0,0), o5s2(20000, 2500), 0653 (40000, 5000), 
0654(60000,0), 0655(60000,5000)) which represent a certain area near the airport. 
We obtained an optimal solution, the obtained noise level is 49.65 dB, the accuracy 
of the results is still very high {max{f.o, e.o) = 5.89e — 07) and the algorithm takes 
no more 130s on a standard PC. This trajectory is about 5 dB less that Ji(54.18 
dB). 
The trajectory characteristics are given in the following figures: 
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Figure 6. Trajectory in 3D with several observers 
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Figure 7. Solutions of (OCP)! 



Almost all state and control variables (except a and V) present large constant 
stages. The control variables Sx and fi are bang-bang between their prescribed 
bounds. 
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5. Conclusion 

We have performed a numerical computation of the optimal control issue. An 
optimal solution of the discretized problem is found with a very high accuracy. A 
noise reduction is obtained during the phase approach by considering the config- 
uration of one and several observers. The trajectory obtained presents interesting 
characteristics and performances. 

Extensions of the analysis on the current problem should include other source of 
noise. This feature is particularly important since improved noise model that better 
represent individual noise sources (engines, airframe,...). It should be remembered 
that this model is focused on single event fiight. Additional researches are needed 
to fully assess the influence of wind and other atmospheric conditions on noise pre- 
diction process. The noise studies have, as yet, been Hmited to a single aircraft 
type equipped with two engines. Further research should consider multiple flights 
model considering airport capacity and nearby configuration. 
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